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We consider the iron pnictides in terms of a proximity to a Mott insulator. The superexchange 
interactions contain competing nearest-neighbor and next-nearest-neighbor components. In the 
undoped parent compound, these frustrated interactions lead to a two-sublattice collinear antiferro- 
magnet (each sublattice forming a Neel ordering), with a reduced magnitude for the ordered moment. 
Electron or hole doping, together with the frustration effect, suppresses the magnetic ordering and 
allows a superconducting state. The exchange interactions favor a d-wave superconducting order 
parameter; in the notation appropriate for the Fe square lattice, its orbital symmetry is d xy . A 
number of existing and future experiments are discussed in light of the theoretical considerations. 
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Introduction: High T c superconductivity has recently 
been discovered in the iron pnictides, with the F-dopcd 
LaOFcAs being the prototype [lj]. Variations include P 
replacement for As, Ni replacement for Fe, and rare-earth 
replacements for La 0, y|, 0, H, @| ■ In addition, F sub- 
stitution for O, which adds itinerant electron carriers to 
the system, could be replaced by, e.c., Sr substitution for 
La, which introduces hole doping [7[. 

Like the cuprates, the iron pnictides have a layered 
structure. The FeAs unit appears to contain all the elec- 
tronic states near the Fermi energy, similar to the case 
of the CuC>2 layer in the cuprates. The electrons par- 
tially occupying the d-orbitals of the iron sites can be 
strongly correlated, as are those on the copper sites in 
the cuprates. At the same time, there are also impor- 
tant differences between the two classes of materials. In 
this letter, we consider the consequences of the unique as- 
pects of the electronic states of the iron pnictides. We will 
frame our discussion in terms of the F-doped LaOFeAs 
family, and touch upon their cousin compounds where 
appropriate. 

One basic question is whether the Mott insulating 
physics plays any significant role in LaOFeAs. The an- 
swer is not necessarily affirmative; there are, for in- 
stance, indications from bandstructure calculations that 
covalency is sizable not only in LaOFeP Q but also in 
LaOFeAs Q. Nonetheless, we argue that there are in- 
direct evidences for the case that LaOFeAs is in prox- 
imity to a Mott insulator. First, the measured electrical 
resistivity is very large, p ~ 5mf2cm at room temper- 
ature [J]. This corresponds to a normalized mean free 
path kpl w hc/e 2 p « 0.5 (where c w 8.7A is the lat- 
tice constant along the normal to the FeAs plane, and 
h/e 2 w 26kf2 is the quantum resistance), which quali- 
fies the system as a bad metal. Second, LaOFeP, which 
has a smaller lattice constant (c = 8.5 A), thus a larger 
internal pressure, is expected to have a larger effective 
bandwidth W but a similar effective Coulomb interaction 



U compared to LaOFeAs. That LaOFeP has a smaller 
U/W compared to LaOFeAs is consistent with the ob- 
servation that LaOFeP is a better metal; its kpi w 1 at 
the room temperature and, indeed, it is superconducting 
with T c w 4 K [2]. These considerations are illustrated 
in Fig. [1] Additional evidence along this direction is pro- 
vided by the lack of a Drude peak in both the measured 
optical conductivity of LaOFeAs , as well as the cal- 



culated one by the DMFT+DFT method [ll|. Thus, our 
approach to LaOFeAs is motivated by important exper- 
imental observations: the very large resistivity and the 
absence of a Drude peak. These are unlike a Fermi liquid 
and therefore are not likely to be accounted for without 
including correlation effects. They imply that most of 
the electronic excitations lie in the incoherent part of the 
spectrum. Even if LaOFeAs is not fully Mott insulating, 
it shouldn't be far away from it. We therefore find it in- 
structive to discuss these systems from a strong-coupling 
point of view. 
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FIG. 1: Placing LaOFeAs and LaOFeP in terms of the control 
parameter U/W, where U is the on-site Coulomb interaction 
and W the effective bandwidth. 

Why is the magnetism so weak? Within this strong 
coupling framework, the issue arises as to why the mag- 
netism is so weak, both in the undoped and lightly-doped 
iron pnictides. Consider first the undoped parent com- 
pound. Valence counting in LaOFeAs yields Fe 2+ , which 
contains six outermost-shell electrons partially filling the 
five 3d orbitals. The degeneracy of the latter is split by 
the crystal field. 
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One characteristic feature seen in the ab initio elec- 
tronic structure calculations is that the splitting among 
the five 3d orbitals is relatively small. Ref. [l2| shows 
that the individual separations among the d levels is on 
the order of, or less than, 0.1 eV. Taking into account the 
typical Coulomb interactions U, of the order of 4-5 eV 
and the Hund's coupling Jh, of the order of 0.7 eV 



we expect the six outermost-shell electrons to occupy the 
3d orbitals in the scheme depicted in Fig. 2a. The asso- 
ciated Mott insulator has spin S = 2. The spin-orbit 
coupling is expected to be considerably weaker than Jh, 
so we shall focus on the spin magnetism. Even if the 
separations between the crystal levels were larger than 
Jh (but still smaller than U), there will still be a double 
degeneracy [HI [3], leading in our consideration to an 
5 = 1 Mott insulator. 

Such a large-spin Mott insulator is expected to be 
strongly magnetic. Yet, a neutron scattering experiment 
[Til has shown that LaOFeAs is an antiferromagnet with 
a rather small ordered moment, on the order of 0.4/is /Fe. 
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FIG. 2: (a) Spin-2 states relevant for the undoped LaOFeAs. 
The crystal levels are according to Ref. The x and y 

we use differ from the standard notation, adopted there, by 
a rotation of 45°; see the main text, (b) Spin-3/2 states that 
become important when electron doping is introduced into 
the FeAs layer. Hole doping will lead to the analogous spin- 
3/2 states, corresponding to five electrons residing on the 3d 
orbitals. 

The issue is even more acute in the doped cases. Elec- 
tron doping will introduce additional states with spin 
3/2, as illustrated in Fig. 2b. Since the relevant local 
states are all magnetic, it is surprising that a relatively 
small amount of doping (say 10% F-doping for O) does 
not preserve the magnetic ordering. 

Magnetic frustration: We propose that the answer 
to these questions lies primarily in magnetic frustration. 
The key feature here is that in the FeAs unit, not only 
do the Fe atoms form a square lattice, but each As atom, 
away from the Fe plane, lies an equal distance from each 
of the four adjacent Fe atoms. Because our focus will be 
on the Fe plane, we find it convenient to use the sym- 
metry classification appropriate for the Fe square lattice: 



we choose the x and y axes to be along the Fe-Fe bond 
direction; these are rotated by 45° from the notation used 
in recent papers. 

Consider the superexchange interactions between the 
3d x 2_ y 2 orbitals of nearby Fe sites. Inspection of the or- 
bitals suggest that the strongest channel of hybridization 
will be with the As Ap x - y or Ap x+y orbital. 

For a pair of ncxt-ncarcst-ncighbor (n.n.n.) Fe 3d x 2_ y 2 
spins, the lowest-energy intermediate state mediating the 
superexchange interaction corresponds to two electrons 
occupying the same Ap x - y (or Ap x+y ) orbital (Fig. 3b). 
Since this intermediate state is a singlet state, the result- 
ing exchange interaction is antiferromagnetic: 
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where V x 2_ y 2 is the hybridization matrix between the Fe 
3d x 2_ y 2 and As Ap x - y orbitals. 

For a pair of nearest-neighbor ( n.n.) 3d x 2 _ y 2 Fe spins, 
the lowest-energy intermediate states correspond instead 
two electrons occupying a pair of distinct Ap x+y and 
^Px-y orbitals (Fig. 3c). The resulting exchange inter- 
action is ferromagnetic: 



J l « -2V x %_ y2 
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Here Jh, p < is the Hund's coupling between the As 
Ap x - y and Ap x + y orbitals, which favors the triplet inter- 
mediate state over the singlet one. Notice that Jh, p is 
relatively small, we expect that the ferromagnetic term 
is small compared to the antiferromagnetic one. 

The dominating matrix elements of the hybridization 
matrix have been given in Ref. [l2T |. The other rele- 
vant hybridizations involve the 3d x / z and 3d y / z orbitals. 
(Here prime denotes the crystallographic axes, which are 
rotated by 45° from the Fe-square-lattice axes.) For 
the n.n.n. interaction, the dominating antiferromagnetic 
terms (with appropriate replacements of the hybridiza- 
tion matrix elements in Eq. ([1])) appear in the diago- 
nal matrix elements and y'z — y'z, as well as 
in the off-diagonal matrix elements x'(y')z — (x 2 — y 2 ). 
For the n.n. exchange J±, the dominating antiferromag- 
netic terms also appear in the diagonal matrix elements 
x'z — x'z and y'z — y'z and in the off-diagonal matrix 
elements x'(y')z — (x 2 — y 2 ). Finally, the n.n.n. exchange 
J2 involves virtual processes associated with only one As 
atom, while the n.n. exchange J\ picks up contributions 
from two As atoms. Taken together, we expect that the 
largest eigenvalues of both the J2 matrix and J\ ma- 
trix correspond to antiferromagnetic interactions, with 
the former larger than half of the latter. 
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dered moment that is considerably smaller than the 
atomic value of order 2/ib- The important point here 
is that it introduces J%l J\ as a tuning parameter, which 
allows for a reduction of the ordered moment. Indeed, as 
J2/J1 is decreased towards a critical value, the ordered 
moment is reduced to zero. The experimental value of 
the moment, as already mentioned, is about 0.4/is/Fe in 



LaOFeAs (15 1 



According to Ref . [l2(] , the 3d — 4p x hybridization ma- 
trix elements and the corresponding energy level separa- 
tions are of the order of 0.8 eV and 1.3 eV, respectively. 
The perturbative expression, Eq. (1), leads to a n.n.n. 
exchange coupling J 2 of the order of 0.5 cV. While it 
is not expected to be quantitatively accurate, the result 
does suggest that the exchange interaction will be sizable. 

For the doped case, the effective model is a matrix 
t — Ji — J2 Hamiltonian: 



Htj — H t + Hj. 
The kinetic component of the Hamiltonian is 



(4) 
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FIG. 3: (Color online) (a) The dominant superexchange inter- 
actions are those between both a pair of nearest-neighbor Fe 
spins (Ji) and a pair of next-nearest-neighbor Fe spins (J2); 

(b) The process contributing to the n.n.n. superexchange 
interaction between a pair of 3d x 2_ y 2 electrons. It involves 
the same <ip x - y orbital, leading to an antiferromagnetic J2; 

(c) The process contributing to the n.n. superexchange in- 
teraction between a pair of 3d x 2 _ y i electrons. It involves 
two orthogonal 4p orbitals (green and magenta, respectively) , 
leading to a ferromagnetic J\. 



The result is the following general form for the spin 
Hamiltonian. 



Hi = 



Jh 



E f 



*i,0, 



(3) 



with mixed J™^ = J" but antiferromagnetic J n .n.n. = 
J% . Here J\ and J2 are both matrices in the orbital ba- 
sis, with matrix elements labeled by a{3. Again, whether 
the local states are spin 2 or spin 1, corresponding to a 
or ft = 1, 2, 3, 4 or a or (3 = 1, 2, depend on whether the 
Hund's coupling Jh is large or small compared to the 
crystal level splittings. 

Eq. ^ specifies a frustrated spin system. In the J2 > 
Ji |/2 case here, the ground state is expected to be a two- 
sublattice collinear antiferromagnet (with each sublatticc 
itself forming a Neel ordering) (lg,ll7|. This spin pattern 
was first proposed for LaOFeAs based on a consideration 
of the Fermi-surface nesting within a spin-density-wave 
picture (HI 18, Tj^. It has subsequently been shown to be 
consistent with the elastic neutron scattering experiment 
in LaOFeAs 

Frustration effects are also important to yield an or- 
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Here, the c Qi j describe constrained fcrmions, which con- 
nect the spin 2 and spin 3/2 configurations at the site i, 
while f^fn = and t n . n . n , — arc the n.n. and n.n.n. 
hybridization matrices. (Recall that a, (3 refer to the d 
orbitals.) The net result of H t is to introduce transitions 
between the spin-2 and spin-3/2 states of the n.n. and 
n.n.n. Fe sites. 

Because the frustration in the superexchange interac- 
tions has already reduced the ordered moment in the un- 
doped parent compound, the magnetic ordering can be 
readily suppressed in the doped materials. This further 
suppression occurs because the F-doping for O induces 
a spin-3/2 (or spin-1/2) substitution of the spin-2 (or 
spin-1) states. Experimentally, the absence of magnetic 
ordering has been shown in LaOi_ x F x FeAs, with electron 
doping of xw 8% 

Superconductivity: Frustration effects, while sup- 
pressing the magnetic ordering, accumulate entropy at 
low temperatures. The relief of this entropy can take 
the form of creating a superconducting order. Precisely 
how this happens is one of the challenging questions in 
strongly correlated electron systems. Still, there are some 
general considerations we can make on the superconduc- 
tivity. 

The proximity to a Mott insulator disfavors isotropic 
s-wave order parameter for the superconducting state. 
Given that the n.n.n. antiferromagnetic exchange in- 
teraction plays a dominant role in the magnetic order- 
ing at half-filling, it is natural that the superconducting 
state has a d xy orbital symmetry. (A mixed pairing state 



(d x 2 



id xy ) may also appear for certain range of J2/J1 



4 



[211].) To see this explicitly, we carry through a simpli- 
fied analysis on a square plaquette. Our consideration 
parallels that of Ref. Q for the cuprate case. 



FIG. 4: A square plaquette showing a spin arrangement of 
the two-sublattice collinear antiferromagnet. 

The square plaquette we consider is illustrated in 
Fig. [4] We will start from the spin-1/2 case. The four- 
electron (N = 4) ground state is well approximated by 



|JV = 4, ps )c<(4 )T c| i T c 2a c u + ' 5 - ii -)|o> 



(6) 



where |0) is the vacuum state, and S.R. denotes spin- 
reversal. 

The N = 2 ground state, at the same time, is well ap- 
proximated by a Gutzwiller-projcctcd Slater determinant 
of two electrons, which can be written in real space as 



interesting subjects. Such studies must await more ac- 
curate determinations from ab initio calculations and/or 
experiment of all the underlying matrix elements [e.g. 
Eqs. (1,2)] and are beyond the scope of this letter. 

A key test for our picture is to experimentally deter- 
mine the relevant spin states in both the parent and 
doped systems, as well as to measure the exchange inter- 
actions from, say, the spin-wave spectra in the undoped 
parent compounds. Studying additional families of mate- 
rials which, in the undoped case, can be placed along the 
U/W axis (Fig. [T]), will allow a fuller exploration of the 
half-filled phase diagram and its doped counterpart. This 
is especially important for the parts of the phase diagram 
that are either deep inside the Mott insulating phase, 
or well into the large-fc^f non-superconducting metallic 
regime. 

Immediately before this paper was finalized, we learnt 
of the work of T. Yildirim [231 ] . which independently con- 
sidered the frustration effect using ab initio calculations, 
and which, in contrast to many other density functional- 
based calculations, gave an ordered moment for the par- 
ent compound in the experimental range. While we be- 
lieve that a consideration of correlation effects is essen- 
tial to account for the "bad metal" properties described 
in the Introduction, there could be other routes to the 
existence of a small ordered moment itself. 

We thank E. Morosan, C. Broholm, H.-P. Cheng, P. 
Dai, K. Haulc, P. Hirschfcld, G. Kotliar for useful discus- 
sions, and acknowledge the support partially provided 
by the NSF Grant No. DMR-0706625 and the Robert A. 
Welch Foundation. 
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where Pq eliminates any double occupancy of a site. It 
is straightforward to show that the pairing operator that 
has the maximum matrix element between \N = 4,gs) 
and \N = 2,<?s) is the d xy pairing operator, 



^(sin k x sin /c. y )c k ,TC-k,] 



(8) 



The construction of the exact \N = 4, gs) and \N — 2, gs) 
for the spin=l/2 well as the equivalent calcula- 

tions for the higher spin cases, can be readily done numer- 
ically. The above suggests that the magnetic exchange 
interactions will promote the d xy pairing, regardless of 
the specific mechanism with which the exchange interac- 
tions cause superconductivity. 

Conclusion: A complete analysis of the matrix 
t — J\ — J i model that we have introduced could reveal 
the existence of other competing phases, such as inho- 
mogencous magnetic structures. The description of the 
phase diagram and its evolution with doping are thus 
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